Experimental observation of shear thickening oscillation 
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We report experimental observation of the shear thickening oscillation, i.e. the spontaneous 
macroscopic oscillation in the shear flow of severe shear thickening fluid. The shear thickening 
oscillation is caused by the interplay between the fluid dynamics and the shear thickening, and has 
been predicted theoretically by the present authors using a phenomenological fluid dynamics model 
for the dilatant fluid, but never been reported experimentally. Using a density-matched starch- water 
mixture, in the cylindrical shear flow of a few centimeters flow width, we observed strong vibrations 
of the frequency around 20 Hz, which is consistent with our theoretical prediction. 



PACS numbers: 83.80.Hj,83.60.Rs, 83.10.Ff,83.60.Wc 

Introduction: Fluid media such as dense colloids or 
dense granule-fluid mixtures are often called dilatant flu- 
ids and known to show severe shear thickening, i.e. its 
viscosity changes discontinuously by orders of magnitude 
upon increasing the shear rate[i|-[3. This is a source of a 
variety of unintuitive curious behaviors of the media [5|- 
0], and might be used for interesting application like a 
body armor 0, but often causes problems in industrial 
situations, such as the damage of mixer motors due to 
over loading [Ij. 

Although the severe shear thickening is a property that 
can be demonstrated easily with common material like 
starch and water, physicists have not reached a general 
agreement on its microscopic mechanism. Originally, the 
order-disorder transition of the dispersed particles was 
proposed [ll lol-llH, but later the hydrodynamics cluster 
formation by the lubrication force is considered to be 
more consistent with experimental observations and nu- 
merical simulations |12l-[l5| . More recently, the transition 
is discussed in connection with the jamming and/or the 
compaction of granules[16l-[T9|. 

Recently, the present authors constructed a fluid dy- 
namics model for the dilatant fluid and found that, in 
a certain shear flow situations, the medium should show 
a spontaneous oscillation, which we call the shear thick- 
ening oscillation^^^^^. Our model is based upon fluid 
dynamics with a phenomenological field (j){r) represent- 
ing an internal state of the medium. The internal state 
(j){r) is determined by the local stress, and the viscosity 
of the medium is a function of ^(r). By assuming sim- 
ple functional relations among these quantities, we have 
demonstrated that the model shows basic properties of 
the dilatant fluids, such as a discontinuous increase and 
hysteresis in the viscosity upon changing the shear rate. 

The shear thickening oscillation that the model pre- 
dicts for the fluid in the shear flow is a salient property, 
but has not been reported in the literature. When the ex- 
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FIG. 1: Schematic illustration of the experimental setup con- 
sisting of a cylindrical container, a rotating rod, and a rotary 
encoder. 



ternal stress is in a certain range, the internal state of the 
medium should alternate between a thickened state and a 
relaxed thin state due to the interplay between the fluid 
dynamics and the internal state dynamics. Such a hy- 
drodynamics oscillation of the medium has neither been 
predicted theoretically nor been reported experimentally. 
By estimating the parameters in the model using exper- 
imental data for the cornstarch- water mixture [2|, the os- 
cillation frequency is expected to be of the order of 10 Hz 
when the flow width is several centimeters, thus it is re- 
ally a hydrodynamic oscillation, not due to a discreteness 
of granules. It should be noted that the oscillation is sup- 
pressed by the viscosity when the flow width is of the or- 
der of millimeters, i.e. a typical size of rheometer sample 
for precision measurement; This may be the reason why 
the shear thickening oscillation has not been reported 
yet. In this paper, we report experimental observation of 
spontaneous oscillation in the cylindrical Couette flow of 
the starch- water mixture. 

Experimental Setup: Figure [1] illustrates our experi- 
mental setup. The fluid flows in the Taylor- Couette ge- 
ometry between the outer cylinder and the coaxial rod 
at the center. The cylindrical container is 24 cm tall. 
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FIG. 2: (Color online) Time evolutions of angular speed of the rod for the flow of h 
stresses. Suspension densities are (a) 41 wt% and (b) 42.5 wt%. 
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but is filled with the fluid upto 22 cm deep. The di- 
ameter of the center rod is 5 cm, while we use several 
outer cylinders with different inner diameters, so that we 
can have the flow of the width h = 1 ^ 5 cm thick- 
nesses. The outer cylinder and the rod are both acrylic, 
and their surfaces are lined with water-proof sand paper 
in order to enforce the no-slip boundary condition. The 
center rod rotates under an externally applied constant 
torque with the outer cylinder being flxed. The exter- 
nal torque is applied by a weight through a steel wire 
wound on the rod; the weight is hung through a spring 
and dumper system (Samini Co., Ltd.), and its weight 
is in the range of 0.5 ~ 10 kg, which gives the external 
stress = 0.14 ~ 2.8 kPa at the surface of the rod. 
The angular speed of the center rod u is recorded by a 
rotary encoder (SP-405ZA, Ono Sokki). As for the fluid, 
we use the suspension of potatostarch (Hokuren) density- 
matched with the aqueous solution of CsCl of the density 
1.6 g/cm^. The upper surface of the fluid is open to the 
air although the container is capped with an aluminum 
plate. 

Shear Thickening Oscillation: Figure [2] shows the 
time evolutions of the angular speed of the rod for the 
potatostarch suspension of the concentrations 41 wt% (a) 
and 42.5 wt% (b) with the flow width h = 4 cm. The 
plots show the oscillations for the various external shear 
stresses Se at the surface of the rod. For both of the con- 
centrations, clear oscillations about 20 Hz are observed 
for the external stress Se ^ 0.2 kPa. The shape of the 
oscillations for the 41 wt% suspension is rather symmet- 
ric and sinusoidal, while for the 42.5 wt% suspension the 
oscillation is nearly symmetric at = 0.2 kPa, but be- 
comes more asymmetric for larger /Se, i.e. a sawtooth 
oscillation that consists of a gradual increase of the an- 
gular speed followed by a sudden drop. For the small 
external stress of = 0.1 kPa, the 20 Hz oscillation dis- 
appears, but a slower oscillation around 5 Hz is observed 
in the case of the thinner fluid of the 41 wt% concentra- 
tion (Fig.2(a)). For the large external stress, we expect 
that the rod will be stuck due to the severe shear thick- 
ening, but in the present experiment, for Se ^ 2.0 kPa, 
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FIG. 3: (Color online) The frequency of the shear thicken- 
ing oscillation for (a) 41 wt% and (b) 42.5 wt% potatostarch 
suspension as a function of flow thickness h. The plots are 
evaluated from the average interval of neighboring peaks, and 
the errorbars represent standard deviation. 



the rod starts slipping suddenly after initial transient. 

Figure [3] shows the oscillation frequency as a function 
of the flow width h for various external stress Se- One 
can see that the frequencies are around 20 Hz for all the 
cases, and we could not flnd any systematic dependence 
on neither h nor Se- For Se higher than 1.0 kPa, the 
center rod starts slipping and we could not obtain clean 
data, but we flnd no clear sign of a systematic change in 
the frequency upto Se = 2.3 kPa. 

Figure [H shows the maximum and minimum angular 
speed of the rod, c^max and cjmin, during a single cycle of 
oscillation against for some values of the flow width 
h for both of the concentrations. Each data point rep- 
resents an average over a single run of the experiment, 
which contains typically 10^ cycles of the oscillation. It 
should be notable that both c^max and cjmin barely de- 
pend on the flow width h. As for the maximum angular 
speed cjmax, the plots show almost linear dependence on 
the external shear stress 5'ei22i]. On the other hand, the 
minimum angular speed c^min stays roughly constant for 
Se ^ 0.5 kPa after the initial decrease for smaller Se- In 
the case of the 42.5 wt% suspension, the oscillation disap- 
pears around Se ^ 0.2 kPa, for which cases the terminal 
angular speeds are plotted with open marks in Fig. 4(c) 
and (d). 
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FIG. 4: (Color online) The maximum and minimum angular 
speed of the center rod during the oscillations as a function 
of the applied stress Se for (a,b) 41 wt% and (c,d) 42.5 wt 
% potatostarch suspensions. Each data point represents an 
average over a single run of experiment, which contains typi- 
cally 10^ cycles of the oscillation. The open symbols in (c,d) 
represent the angular speed of the non-oscillating flows for 
small Se- 



The linear dependence of cjmax on Se leads us to intro- 
duce an average viscosity 77reiax for the flow during the 
relaxed state. If we assume the steady cylindrical Cou- 
ette flow with a constant angular velocity cj, the viscosity 
would be given by 



(1) 



where Ri = 2.5 cm is the radius of the center rod and 
R2 = Ri -\- h is the radius of the outer cylinder. For 
the 42.5wt% suspension in the relaxed state, the plots 
in Fig. 4(c) shows c^max ~ IS^e rad/s. Substituting this 
^max into u in Eq.(l), we have r/reiax = 20 ~ 30 Pa-s, 
which is consistent with the previous viscosity measure- 
ment before the thickening transition [2]. For the 41 wt% 
suspension, Fig|4l^a) shows cjmax ^ lO+ll^e rad/s, which 
has a substantial non-zero extrapolation at = 0. This 
may be due to the inertia effect in the oscillatory state 
and our steady state formula, Eq.([T]), is not valid for this 
case. 

In the thickened state, c^min is nearly constant around 



0.5 rad/s for 0.5 < 5'e < 1.5 kPa for the 42.5wt% sus- 
pension, thus the viscosity in the thickened state may be 
estimated to be of the order of 10^ Pa-s, which seems to 
be somewhat smaller than the value observed in Ref.[2]. 
This may suggest that the whole system is not thickened 
uniformly, but the thickening region is localized as we 
will discuss later. 

Comparison with model behaviors: Some features we 
observed are consistent with those predicted by our phe- 
nomenological model [21]. The shear flow starts oscillat- 
ing under a constant shear stress larger than a certain 
value. The oscillation is sinusoidal for the stress near 
the threshold, but becomes asymmetric for larger stress 
in the sawtooth shape. The observed frequency of the 
oscillation is about 20 Hz, which is roughly in the range 
expected by the model for the flow width of a several 
centimeters. 

On the other hand, there are some other features that 
seem difficult to interpret in terms of the simple shear 
flow oscillation. The oscillation frequency does not de- 
pend significantly neither on the external stress Se nor 
on the flow width /i, whereas the simulations of cylindri- 
cally symmetric oscillation show that the frequency tends 
to become slightly higher as Se increases, and that the 
frequency decreases signiflcantly as the flow width h in- 
creases. It is also tricky to interpret the behavior of the 
maximum cjmax and minimum angular speed cjmin of the 
oscillation: cjmax is proportional to Se while cjmin stays 
almost constant, and neither of them do not depend on 
the flow width h. We flnd that these features are not easy 
to reproduce by the model as long as we assume that the 
flow is cylindrically symmetric. 

In fact, the simulations in two dimensions show that a 
localized thickened band region appears when the cylin- 
drical Poiseuille flow is destabilized beyond a certain Se- 
Figure [5] shows the results of numerical simulations of the 
fluid dynamics model for the time developments of 

the flow when the inner cylinder is driven toward CCW 
with the outer cylinder being flxed. The velocity fleld 
v{r) and the spatial variation of the viscosity r]{r) are 
shown by the arrows and the gray scale maps, respec- 
tively. 

The initial destabilization process of the cylindrically 
symmetric flow is shown in Fig|5ja) after the system 
starts from the static state. As the inner cylinder starts 
rotating, a uniform circular flow is induced and the region 
near the inner cylinder thickens, then a cylindrically sym- 
metric state is destabilized. A localized thickened region 
forms a band structure to extend in the radial direction, 
then the flow slows down drastically when the thickening 
band reaches the outer cylinder. 

After the initial transient, the system undergoes steady 
periodic oscillation, where the thickening band breaks off 
and extends periodically during the course of the oscil- 
lation (Fig|5jb)). The band goes around with the inner 
cylinder. Our tentative results of the simulations show 
that, for this oscillatory flow with a thickening band, 
the frequency is determined by the way how the band 
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FIG. 5: Time developments of cylindrical shear flows by the simulations of the fluid dynamics model of shear thickening pH: 
(a) the initial destabilization of the cylindrically symmetric flow from the static state, (b) one cycle of steady oscillation with 
a localized thickening band. The inner cylinder is driven toward CCW by the external shear stress Se with the outer cylinder 
being fixed. The rotation of the inner cylinder is represented by the solid lines at the center. The velocity field v{r) and the 
spatial variation of the viscosity 77 (r) relative to that in the relaxed state 770 are shown by the arrows and the grey scale maps, 
respectively. The simulation parameters are the driving shear stress 5*6 = 3, the radius of inner cylinder i?i = 1, and the flow 
width /i = 2 in the dimensionless unit system defined in Ref.^]. The model parameters are r = 0.1, A = 1 and (^m = 0.85. 
The moment of inertia of the inner cylinder is neglected. 



breaks off, and its dependence on the flow width h and 
the shear stress Se is weak in comparison with the case of 
the cylindrically symmetric flow. This may explain our 
experimental data. Detailed behavior of the oscillatory 
flow with this localized thickening band is now under in- 
vestigation. 

In summary, we observed the shear thickening oscilla- 
tion around 20 Hz in the shear flow of the Taylor- Couette 
geometry with several centimeter flow width as has been 
predicted theoretically. The frequency and the amplitude 



of the oscillation barely depend on the flow thickness. 
Numerical simulations in two dimensions suggests that 
the oscillation occurs with a thickening band that breaks 
off and extends periodically. 
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